Non-Abelian spin-singlet quantum Hall states: 
wave functions and quasihole state counting 
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We investigate a class of non-Abelian spin-singlet (NASS) quantum Hall phases, proposed pre- 
viously. The trial ground and quasihole excited states are exact eigenstates of certain k + 1-body 
interaction Hamiltonians. The k = 1 cases are the familiar Halperin Abelian spin-singlet states. 
We present closed-form expressions for the many-body wave functions of the ground states, which 
for k > 1 were previously defined only in terms of correlators in specific conformal field theories. 
The states contain clusters of k electrons, each cluster having either all spins up, or all spins down. 
The ground states are non-degenerate, while the quasihole excitations over these states show char- 
acteristic degeneracies, which give rise to non-Abelian braid statistics. Using conformal field theory 
methods, we derive counting rules that determine the degeneracies in a spherical geometry. The 
results are checked against explicit numerical diagonalization studies for small numbers of particles 
on the sphere. 

PACS: 73.43.-f, 71.10.-w, 71.10.Pm 
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I. INTRODUCTION AND SUMMARY 

The observation jl]-|| of a quantum Hall (QH) state 
at an even- denominator filling factor, v = 5/2, stimu- 
lated the development of trial wave functions outside the 
usual hierarchy (or later, composite fcrmion) approach, 
which generates only odd-denominator fractions. The 
5/2 state is interpreted as half- filling of the first excited 
Landau level (LL), the lowest one being filled with elec- 
trons of both spins, and can be mapped to half-filling of 
the lowest LL, with a suitable Hamiltonian. There are 
now strong indications that this state is spin-polarized 
P,|J, and described by the paired "Pfaffian" state of 
Moore and Read (MR) g], which has filling factor 1/2. 
This state was originally proposed as an example that 
manifests non-Abelian braid statistics of its quasiparticle 
excitations [||. Generalizations exist in which the parti- 
cles are "clustered" in fc-plets (k — 1, 2, 3, ... ), but 
still spin-polarized 0. In these states, the non-Abelian 
statistics are associated with parafermion conformal field 
theories (CFTs). 

It is well known that, despite the presence of a strong 
magnetic field, spin-singlet QH states are sometimes fa- 
vored over their spin-polarized counterparts. The possi- 
bility to manipulate the Zeeman splitting by the appli- 
cation hydrostatic pressure and by tilting the magnetic 
field has opened the possibility of systematic studies of 
transitions between spin-polarized and non-polarized QH 
states (see, e.g., Q). In this light, states analogous to 
the spin-polarized clustered QH states of j?J , but with a 
spin-singlet structure, have been constructed M. In ||, 
trial wave functions for these non-Abelian spin-singlet 
(NASS) states have been written in terms of correlators 
in a CFT describing parafermions associated to the Lie 



algebra SU(3). We remark that one may consider an al- 
ternative series of NASS states, whose algebraic structure 
is related to SO (5) rather than SU(3). The simplest state 
of this type is a paired spin-singlet state that exhibits a 
separation of spin and charge in the quasihole excitation 
spectrum The SO(5)-based NASS states will not be 
discussed in this paper. 

In the present paper, we study in detail some of the 
properties of the NASS states, paying special attention 
to the case k = 2. We give explicit closed form expres- 
sions for the ground state wave functions, and study the 
degeneracies of their quasihole excitations. The degen- 
eracies of states with fixed spins and fixed well-separated 
positions of the quasiparticles are the origin of the non- 
Abelian braid statistics. 

We first strengthen the case for the existence of the 
incompressible phases of matter with the universal prop- 
erties of the states of Ardonne and Schoutens (AS) M, 
by showing that their trial wave functions for the ground 
state and for states with quasiholes are exact zero-energy 
eigenstates of certain k + 1-body interaction Hamiltoni- 
ans for particles in a single LL, in a similar way as the 
spin-polarized cases 0. The explicit closed-form wave 
functions for the ground states are obtained. In the 
study of the quasihole degeneracies, we then follow two 
complementary approaches. The first is an analytical 
path, which relies heavily on the formal structure of the 
associated parafermion CFT, and on the analogy with 
earlier studies for spin-polarized non-Abelian QH states 
PI O. While at present we lack explicit expressions 
for the many-body wave functions describing the quasi- 
holes, we have enough control to derive explicit counting 
formulas for the degeneracies, for k = 2, for particles 
on a sphere. The second approach is a numerical study 
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of the k + 1-body Hamiltonian for the case k — 2, on 
the sphere. The numbers of zero-energy states for each 
number of electrons N and of quasiholes n considered 
are in exact agreement with the analytical derivation. In 
addition, we study the excitation spectrum of the same 
Hamiltonian, and compare the ground state with that of 
electrons interacting via the lowest LL Coulomb interac- 
tion. 

A highlight of the analytical approach in this paper 
comes in the derivation of the total degeneracies of quasi- 
hole states. In the CFT set up (which will be described in 
more detail in Section [n]) the QH states are described as 
conformal blocks of "particle" and "quasihole" operators. 
The particle operator factorizes as a product of a vertex 
operator and a parafermion field, and the quasihole op- 
erator is the product of a vertex operator and a so-called 
spin field of the parafermion CFT. The non-trivial fusion 
rules of these spin fields cause a degeneracy of the ground 
states in the presence of quasiholes at fixed positions and 
spins. There is also further degeneracy associated with 
the positions and spins of the quasiholes, which is finite 
on the sphere. The two contributions need to be com- 
bined in the right way. It turns out, as in earlier cases, 
that the various states which stem from the non-trivial 
fusion rules have a different spatial degeneracy. There- 
fore, we can not just multiply the two degeneracy factors, 
but we need to break up the degeneracies due to the fu- 
sion rules. To accomplish this task, and arrive at the final 
counting formula, we analyze "truncated chiral spectra" 
in the SU(3) parafermion CFT, using the methods of Q . 

This paper is organized as follows. In Section p], we 
explain in which way CFT is used to describe QH states, 
and review the NASS states as correlators. In Section 
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we introduce the k + 1-body interaction, and show 
that the AS correlators give zero-energy ground states. 
In Section |[v]we give the explicit ground state wave func- 
tions for the NASS states, and discuss their spin-singlet 
properties. Section [v| describes the correlators which give 
the states with quasiholes present. The derivation of the 
counting formula is done in Sections 
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t o [1X| , using the 

method which is outlined above, with Eq. ( |49|) as the final 
outcome. Explicit results of this formula for the degen- 
eracy of the ground states in the presence of quasiholes 
are given in the same Section for several N (the number 
of electrons) and n (the number of quasiholes). In Sec- 
tion [x] we present numerical diagonalization studies on 
a sphere, finding full agreement with the analytical ex- 
pression obtained in Section IX, and compare the states 



with the ground state of the Coulomb interaction. 



II. QHE-CFT CORRESPONDENCE 

In the QHE, following Laughlin p5[ , trial wave func- 
tions have long been used as paradigms that represent an 
entire phase of incompressible behavior. This notion was 



reviewed in Ref . |t]] , so in this paper we will concentrate 
on the properties of trial states and their position-space 
wave functions. As explained by MR fj], many QH trial 
wave functions of the 2 + 1 dimensional system can be 
obtained as conformal blocks (i.e., chiral parts of correla- 
tion functions or "correlators" ) in a suitable chiral CFT 
in 2 Euclidean dimensions, as we will briefly explain. 

For particles with complex coordinates Zj = Xj + iyj , 
j = 1, ... , for N particles, we will use reduced 
wave functions ^(z), and neglect spin temporarily. For 
the lowest Landau level (LLL), the reduced wave func- 
tion must be holomorphic in the z^'s. For particles in 
the plane, the full wave function ^(z) is recovered by 
multiplying ^(z) by exp(— J^i l z i| 2 /4); we have set the 
magnetic length to 1. For particles on the sphere |l6[ , 
the coordinate zi refers to stereographic projection, and 
the full wave function is recovered by multiplying by 
Uii 1 + \z l \ 2 /^R 2 )~ (1+N '* ,2 \ where is the total num- 
ber of magnetic flux quanta through the sphere Jl2],[7[. 
In the latter case, the reduced wave function ^(z) must 
be a polynomial of degree no higher than N$ in each z^, 
so that the z component of angular momentum of each 
particle lies between and —N^/2. Note that we use 

the term "particles" for the underlying particles, which 
are either charged bosons or charged fermions (electrons), 
because it is convenient to consider cases of either statis- 
tics together. 

The simplest example of a state with a uniform density 
(a state of zero total angular momentum on the sphere 
fll6f ) is the Laughlin wave function [fl5| : 
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M 



(1) 



for a fixed integer M . The filling factor can be defined 
for a sequence of states as v = limjv^oo N/N^ 7 where 
is identified with the largest power of any zi in the state. 
For the Laughlin state it is v = -h. Note that this func- 
tion is antisymmetric (describes fermions) for M odd, 
and symmetric (describes bosons) for M even. 
The Laughlin wave function can be obtained as 



lim (z oo r{V e (z 1 )...V e {z N )e-^ N ^(z oc )) , (2) 



with V e (z) — exp(ivM</j) a chiral vertex operator in 
the c = 1 chiral CFT describing a single scalar field 
ip compactified on a radius R 2 = M. The operator 
e -iVMNtp^ z ^ brings in a positive background charge, 
which guarantees the overall neutrality of the system. 
The constant a must be chosen in such a way that the 
effect of the background charge does not go to zero in 
the limit z^ — > oo; in Eq. (j^) we need a = MN 2 . This 
procedure is simpler for our purposes than the uniform 
background charge used in MR H , though the latter has 
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the additional feature of reproducing the factors in the 
unreduced wave function. 

Other, more complicated, QH states can be con- 
structed by invoking more complicated CFTs. The CFT 
framework guarantees that a number of consistency re- 
quirements for such states are met Q. The trial wave 
functions become more meaningful, and the correspond- 
ing phase really exists, when there is a (local) Hamilto- 
nian for which the trial state is the nondegenerate ground 
state, and the excitation spectrum (for the same N,p as 
the ground state) has a gap in the thermodynamic limit. 
Short range 2-body interaction Hamiltonians with these 
properties for the Laughlin state were found by Haldane 
ijlrjf , and 3-body interactions for which the MR state is an 
exact zero energy eigenstate were found beginning from 
the work of Ref . jnj . Read and Rezayi (RR) discovered 
that these constructions are the first two cases in an 
infinite sequence, and found the parafcrmion states as the 
exact zero energy eigenstates of k + 1-body interactions 
for all k. Here we will show similarly that, in the case of 
particles with spin, the NASS states of Ref. || are exact 
eigenstates of zero energy for k+l body spin-independent 
interactions, with the Halperin state and the correspond- 
ing 2-body interaction |l6| as the only case known previ- 
ously. First, we recall the construction in Ref. S, then in 
the following Section we establish that the wave function 
defined by a correlator here is a zero-energy eigenstate of 
a k + 1-body interaction Hamiltonian. 

The NASS states proposed in M can be viewed as 
non-Abelian generalizations of the Abelian spin-singlet 
Halperin states labeled as (m + 1, m + 1, m) (see Eq. 
(^) below) , or alternatively, as spin-singlet analogs of the 
spin-polarized "clustered" or "parafermion" states of RR 
0. The filling fraction of the NASS states is v = 2k lj +3 , 
with M an integer (M is odd when the particles are 
fcrmions, even when they are bosons). The wave func- 
tions of these states are constructed as conformal blocks 
in basically the same way as was done above for the 
Laughlin state. In the basic case M — 0, the component 
of the ground state with any set of N/2 of the particles 
having spin up, and the remainder spin down, is defined 
(up to phases that may be needed when reconstructing 
the full state from these components) as the correlator 



k,M=0 
NASS 



3 jy2 TV 

lim (zoo) 2k (exp(i — — (a 2 -a 1 )-0)(z oc> ) 
oo^oo 2Vk 



)> • (3) 

In this equation the "particle operators" (to avoid con- 



fusion, we emphasize that this means the operators that 
correspond to the particles in the CFT, not the opera- 
tors that create actual particles in the 2 + 1-dimensional 
system) are current operators B a (z) in an SU(3)fc (i.e., 
level k) Wess-Zumino-Witten CFT. These currents can 
be written in terms of two bosons (p ~ {tpx^ip^) and 
a Gepner parafermion i/j a associated to SU(3)fe/[U(l)] 2 
p8[ . The currents are labeled by the corresponding roots 
a of SU(3) 



B a (z) = tp a exp(?a • (pj Vk){z) 



(4) 



The roots are given by a.\ = (V2, 0), a 2 — 
(— -\/2/2, "v/6/2). We note that these two roots form an 
SU(2) doublet under an SU(2) subalgebra of SU(3); the 
embedding of the subalgebra is isomorphic to that given 
in terms of 3 x 3 Hermitian matrices [generators of SU(3)] 
by the 2x2 Hermitian blocks at the upper left corner. 
For such an embedding, there is also a U(l) subalgebra 
[generated by "hypercharge" diag(l, 1, —2); the particles 
carry hypercharge 1] that commutes with the SU(2) and 
corresponds to the particle number. Note that a (hyper- 
charge at infinity is again needed for neutrality. 

Working out the vertex-operator part of this correlator, 
we arrive at the following factorized form of the NASS 
state (after multiplication by an additional Laughlin fac- 
tor to obtain general M 1 ) 



T J. 

NASS^l ) • • • ) z N/2' Z l ' ■ ' ■ ' / -N/2> 
(t/j ai (4) ■ ■ ■ OLo^-M (4) ■ ■ ■ ( Z N/ 2 )) 



l) 



"N/2' 
1/k 



K( z J;4 



(5) 



The wave function ^ 
functions lfl9| 



(2,2,1) 



is one of the Halperin wave 



{m',m',m) ( t 

1^1 ! 



T . 1 

' Z Af/2' Z l ' ' 



"N/2> 



(G) 



l<3 



i<j 



1,3 



The latter give rise to spin-singlet states whenever m! — 
m + 1 (l(| . The wave function Eq. (||) contains a term 
which is a correlator of parafcrmions, the explicit form 
of which will be found below. 

We also mention here that the CFT construction im- 
plies that the number of sectors of edge states (represen- 
tations of the chiral algebra), and the number of ground 



1 We can also obtain the Laughlin factor by using the parti- 
cle operators (|l!j) and (^) in the correlator, together with a 
suitably adjusted background charge. 
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states (conformal blocks with N particle operators in- 
serted) on the torus for N divisible by 2k, are both given 
by (k + l)(k + 2)(2fcM + 3)/6, which is an integer. For 
M = 0, this coincides with the numbers for SU(3)fc cur- 
rent algebra. The filling factor is P/Q = 2k/(2kM + 3), 
so if P and Q are defined as being coprime, then the de- 
nominator Q — 2k M + 3 unless 2k and 3 have a common 
factor, that is unless k is divisible by 3, in which case we 
have Q — (2kM + 3)/3. The number of ground states on 
the torus is then always divisible by the denominator Q 
of the filling factor, as it should be. We also expect that 
for some k values there are ground states on the torus for 
other N values, as for the MR states [^0[ and RR states. 

III. SOLUTION OF K + 1-BODY HAMILTONIAN 

It is known |l(| that the Abelian Halperin spin-singlet 
state is the unique (on the sphere) exact zero-energy 
eigenstate of a two-body interaction Hamiltonian. Other 
than the trivial case (1,1,0), which is two filled Landau 
levels, where the Hamiltonian in question is zero, the sim- 
plest case is m = 1 in Eq. (^J), the (2,2,1) state, which 
corresponds to M — in Eq. (||) or (||). In the latter 
case this Hamiltonian is simply a repulsive <5-function in- 
teraction between any two particles. As in Ref. fll2|,|7j], 
it is simplest to start by generalizing this M = case. 
Because higher M values are obtained by multiplying by 
additional Laughlin factors, the Hamiltonians for M = 
can be straightforwardly extended to M > by extend- 
ing the range of the k + 1-body part, and adding 2-body 
interactions as needed, which have the effect of requiring 
zero-energy states to contain the Laughlin factors. We 
will not detail this here, however, see Sec. g below. 

The natural choice of Hamiltonian for M = is to con- 
sider the k + 1-body ^-function as in Ref. 0, but here 
for particles with spin. The Hamiltonian (acting within 
the LLL) is 

H = V ]T 5 2 ( Zll -z l2 ) 

il<i2<---<ik+l 

x6 2 (z, l2 - z l3 )---S 2 (z lk ~z lk+1 ), (7) 

with V > 0. Note that here we have reverted to label- 
ing the particles independently of their spin. For this 
Hamiltonian, a state is a zero-energy eigenstate if it van- 
ishes whenever any k + 1 particles coincide; for this to be 
satisfied for some nontrivial states, the particles must be 
bosons. 

We will now show that the correlator as in Eq. (||) is 
such a zero-energy eigenstate. The argument we give here 
makes direct use of the current algebra satisfied by the 
currents, and also sheds new light on the previous spin- 
less case of RR, where a slightly different argument was 
used. Without loss of generality, we can consider letting 



the first k + 1 particles, of either spin, come to the same 
point, say z = 0, that is z^ 1 , z% 2 , . . . , all — > 0. 

In the standard radial quantization scheme for CFT, we 
can consider the current operators as acting in a Hilbert 
space that is built starting from a highest weight state 
that in the present case is simply the vacuum |0) of radial 
quantization about the origin z = 0. As zp tend to one 
by one, the resulting operator product expansion (ope) 
contains no singular terms. This follows from the stan- 
dard current-algebra ope's of the currents, together with 
the fact that the roots a,\ and — ui do not sum to either 
or another root (for simplicity, let us replace these two 
roots by the natural notation a =\ and J., respectively). 
Indeed the only nonvanishing term as the z°' tend to zero 
sequentially is the operator at that corresponds to the 
state 

B^-xB^-x ■ ■ ■B ak+1 _i\0) (8) 

in the highest weight representation, and we need to show 
that this vanishes for all choices of ax, ■ ■ ■ , erfc+i- Here 
we have used the modes of the currents, 

B a {z) = ^-™- 1 J B a ,„, (9) 

n 

which holds for all generators a of SU(3), not only roots. 
In fact, the commutation relations of the affine Lie al- 
gebra for these modes also imply that -Bf,-i and £?x,-i 
commute, so we need only to prove 

(B i ,_ 1 r(B T ,_ 1 )' £ + 1 - m |0)-0, (10) 

for all m in the range < m < k + 1 . 

Let us begin by choosing m — 0. Then we need to 
show that 

(i? T ,_ 1 ) fc+1 |0)=0. (11) 

But this is simply the pure-current-algebra null-vector 
equation, which first entered the physics literature in 
Refs. pi| , P2[ . Thus this is satisfied in the irreducible, uni- 
tary vacuum highest weight module of the SU(3) affine 
Lie algebra at level k, and there are similar equations for 
all current algebras, and for each "integrable" highest 
weight representation. This already suffices to rederive 
the RR states, which are related to SU(2) current algebra 
@. The RR states are in fact the same as those in Eq. (^) 
above, but with all spins |, so the only root that appears 
lies in an SU(2) subalgebra (the spin singlet property of 
the above states is of course lost when this is done, and 
the charge at infinity should be adjusted). Hence, we 
have very quickly rederived the fact that the RR states 
for M = are the zero-energy states of the k + 1-body 
(5-function Hamiltonian. 

It is straightforward to complete the proof for the 
NASS states. Essentially, we use the SU(2) symmetry 
under which and f° rm a doublet, and notice 
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that the set of states labeled by m forms a highest weight 
multiplet under this algebra, of SU(2) spin k/2. (This is 
clear from the 2 + 1 point of view, where we are looking 
at states of k + 1 spin 1/2 bosons all at the same point.) 
Then since the highest weight vanishes, all the others do 
also, which completes the proof. 

As in RR, a similar argument also establishes that 
quasihole states, written as similar correlators but with 
spin fields (primary fields of the SU(3) current algebra) 
inserted at the quasihole positions JsJ , are exact zero- 
energy eigenstates. These are considered explicitly in 
Sec. [y] below. Similar arguments also imply that the 
zero-energy ground states of H on the torus are given 
by correlators for some number N of the above fields in- 
serted, with the number of such ground states (for N 
divisible by 2k) already given in Sec. ||. 

IV. GROUND STATE WAVE FUNCTION 

Based on the results of the previous Section, we expect 
the structure of the trial wave functions — that is, of the 
chiral correlators (|J) — of the NASS states to be similar to 
that of the RR states, and also to generalize the Halperin 
(2,2,1) state. The RR wave functions were constructed 
by dividing the (same-spin) particles into clusters of k, 
writing down a product of factors for each pair of clus- 
ters, and finally symmetrizing over all ways of dividing 
the particles into clusters. Hence in the case with spin, we 
guess that we should divide the up particles into groups of 
k, the downs into groups of k and then multiply together 
factors that connect up with up, down with down, or up 
with down clusters, and finally ensure that the function 
is of the correct permutational symmetry type to yield a 
spin-singlet state (in particular, it should be symmetric 
in the coordinates of the up particles, and also in those of 
the downs). We expect that the up-up and down-down 
parts of this should closely resemble the RR wave func- 
tions, before the symmetrization; it was shown in Ref. 
that the functions found there vanish when k + 1 par- 
ticles come to the same point, even inside the sum over 
permutations that symmetrizes the final function. These 
considerations guided the following construction. 

Due to the spin-singlet nature of the state, the wave 
function will be non-zero only if the number of spin-up 
and spin-down particles is the same. Furthermore, there 
must be an integer number of clusters, so the total num- 
ber of particles N must be divisible by 2k, and will be 
written as N = 2kp, where p € N. One example was al- 
ready given in , namely the wave function for the case 
k = 2, M = with the number of particles equal to 4 
(i.e., p = 1), 

^NASS V z l > z 2' z l i Z 2J — 

(4 - 4)(4 - 4) + (4 - 4)(4 - 4) ■ (12) 



This is part of the two-dimensional irreducible represen- 
tation of the permutation group on 4 objects, S4, as can 
easily be seen. This is the correct symmetry type to ob- 
tain a spin-singlet state, as we discuss further below. 

We will now describe the different factors that enter 
the NASS wave functions. Because the only effect of M 
being non-zero is to give an overall Laughlin factor, we 
will assume at first that M = 0. First we give the factors 
that involve particles of the same spin, say spin up. They 
are the same as in RR 0. We will divide the particles 
into clusters of k in the simplest way, 

(zl,... ,z\),(z\ +l ,... ,z\ k ),... , (Z( p _^ k+1 , . . . ,Zp k ) , 

(13) 

and the same for the z^'s. (In a more precise treatment, 
we would say that the first N/2 particles are spin up, 
the remainder spin down.) We write down factors that 
connect the ath with the 6th cluster: 

Xa.b = ( Z (a-l)fc+l — Z (&-l)/c+lX Z {a-l)fc+l _ Z (b-l)k+2) 
X ( 2 (a-l)fc+2 _ Z (fc-l)fc+2)( Z [a-l)fc+2 _ Z (b-l)k+s) 

x • • • x (4k - 4 k )(4k - 4b-i)k+i) • ( 14 ) 

For k = 1, we would write xtb = (4 ~ 4) 2 ■ The factors 
that connect up with down spins are simpler: 

Aa,fc — ^(o-l)fc+l Z (b-l)k+l> 
X ( Z (a-l)k+2 ~ Z (fc-l)fc+2) ' ■ ' (4k ~ 4k) ■ 

For k = 1, the factor would be Xab* = (4 ~ 4)- W° 
multiply all these factors for all pairs of clusters, up-up, 
down-down, or up-down: 

flxib f[x z c y fixif (i6) 

a<b c.d e</ 

Notice that for k = 1, we do obtain the Halperin (2,2,1) 
wave function. 

To obtain a spin-singlet state when the spatial function 
is combined with the spin state (which lies in the tensor 
product of N spins 1/2), some symmetry properties must 
be satisfied. For the M — case, the particles are bosons, 
hence the full wave function must be invariant under per- 
mutations of spins and coordinates of any two particles. 
This can be used to obtain the correct form of the func- 
tion from that component in which, say the first N/2 are 
spin up, the rest spin down, as above, so knowledge of 
that component is sufficient. The requirement that the 
full wave function be a spin-singlet can be shown to re- 
duce to the Fock conditions: the component just defined 
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must be symmetric under permutations of the coordi- 
nates of the up particles, and also of the down particles, 
and must also obey the Fock cyclic condition, as given 
in Ref. (modified in an obvious way for the boson 
case) . These three conditions can be shown to imply that 
the spatial wave function is of a definite permutational 
symmetry type (belongs to a certain irreducible repre- 
sentation of the permutation group), that corresponds to 
the Young diagram with two rows of N/2 boxes each. In 
general, given a function of arbitrary symmetry, a Young 
operator can be constructed that projects it onto a mem- 
ber of the correct representation (though the result may 
vanish); this construction generalizes the familiar sym- 
metrization and antisymmetrization operations. For the 
present case, the Young operator is the following op- 
eration, equivalent to summing over the function with 
various permutations of its arguments, and some sign 
changes: First, antisymmetize in z\, £/\r/2+i> then in z 2l 
zn 72+2; • ■ ■ , z N/2 , z N ; then symmetrize in z\, ... , z N/2 ; 
then finally symmetrize in 2^/2+1, ■ ■ ■ , This clearly 
satisfies the first two requirements of Fock, and can be 
proved to satisfy also the cyclic condition. It remains to 
check that it is nonzero, we believe it is. Incidentally, the 
application of the Young operator is the analog of sym- 
metrizing over the down spins in the spatial wave func- 
tion of the permanent state (see e.g. Ref. to which 
it reduces for the case of BCS paired wave functions of 
spin 1/2 bosons (there are similar statements in the more 
familiar case of spin-singlet pairing of spin-1/2 fermions). 
However, based on the example of the Halperin (k = 1) 
case, we also considered the function defined as in Eq. 
(|l6|), and then simply symmetrized over all the ups and 
over all the downs. For the Halperin function [which in 



fact is already symmetric in Eq. (16)], this satisfies the 



cyclic condition, as can be seen using the fact that the 
(1,1,0) state is a Landau level filled with both spins, plus 
the Pauli exclusion principle for fermions. For k = 2, 
3, we verified the cyclic condition numerically for several 
moderate sizes. Hence, we expect that this simpler form 
actually works for all k (as well as for all N divisible by 
k). Apparently, this procedure and the application of the 
Young operator give the same function in the end (up to 
a normalization). 

For M = 0, our wave function is then: 



A-.O 



,\.\ss - Sym J \ 

a<b 



nz\z l 
Xc,d 



IK/ 



(17) 



where Sym stands for the symmetrization over the ups 
and also over the downs. This function is nonzero, as 
may be seen by letting the up coordinates coincide in 
clusters of k each, and also the downs, all clusters at dif- 
ferent locations, and making use of the result in RR |t]] 
that only one term in the symmetrization is nonzero in 
the limit. This term is the Halperin {2k, 2k, k) function 
for 2p particles. To obtain the wave function for general 



M, we multiply by an overall Laughlin factor, fyjf . 

We can give a simple proof that our wave function (for 
M = 0) vanishes if any k + 1 particles, each of either spin, 
come to the same point. This works also for the RR wave 
functions, and is simpler, though less informative, than 
the proof in RR fij. It works term by term, inside the 
sum over permutations in the symmetrizer. Thus, with- 
out loss of generality, we may use the simple clustering 
considered above. We note that on the clock face formed 
by the labels 1, . . . , k within each cluster, there is always 
a factor connecting any two particles at the same posi- 
tion, regardless of their spin. This factor vanishes when 
the particles coincide. Since there are only k distinct 
positions, when k + 1 particles come to the same point, 
the clock positions must coincide in at least two cases, 
so that the wave function vanishes, which completes the 
proof. 

We do not have a direct general proof of the equality of 
these explicit wave functions and the formal expressions 
Eq. (||), but we have performed a number of consistency 
checks. First, the wave functions arc polynomials of the 
correct degree. From Eq. (||), we can infer what the total 
degree should be. The parafermions of the correlator con- 
tribute with (see pi) -l-2fcp-(l-i). The factors of the 
2,2,1 part are 2 • f^\kp(kp- 1) and 1- 1 • [kp) 2 . Adding 
these gives, for M = 0, pk(3p — 2). We need to check 
whether Eq. @ gives the same degree. For the ith up 
particle, the degree of zj in the product of up-up factors 
is 2(p— 1), and in the up-down factors is p. Thus the net 
degree in zj is = 3p — 2 = 3N/2k — 2, or for general 
M, N (j> = 3p + M(N - 1) - 2 = (M + 3/2k)N -2-M. 
This gives the filling factor v = 2k/(2kM + 3) @, which 
reduces to that for the Halperin states for k — 1, and 
also the shift, defined as N$ = N/v — S, which here is 
S = M + 2 on the sphere (for more on the shift, see 
Ref. [Q). Finally, the total degree is N/2 times that in 
zj , namely kp(3p — 2) for M — 0, the same as for the 
correlator. Also, the numerical work described in Sec. |x] 
below confirms that the ground state of the appropriate 
Hamiltonian on the sphere for k = 2, M — 1 at the given 
number of flux does have a unique spin-zero ground state 
at zero energy, so that the correlator and the wave func- 
tion constructed above must coincide. This also implies 
that the wave functions above must be spin singlet. 



V. CORRELATORS CORRESPONDING TO 
STATES WITH QUASIHOLES 

To obtain wave functions for NASS states with quasi- 
holes, one inserts corresponding operators into the cor- 
relator that corresponds to the ground state wave func- 
tion. Here we will give the form of the correlators, using 
standard CFT techniques, but not the wave functions. 
The operators we insert have the form of a spin field 



G 



times a vertex operator, similar to the RR states |j (note 
that the term "spin field" is traditional, and has no rela- 
tion to the SU(2) "spin" symmetry). In the exponent of 
the vertex operators, the fundamental bosons are multi- 
plied by (fundamental) weights of the Lie algebra SU(3): 
zui = (V2/2, y/6/6), w 2 = (0, y/6/3). We take the quasi- 
holcs from the triplet 3 of SU(3). The corresponding 
operators are 

C roi (w T ) = cr roi exp(ira7i • Lp/Vk)(w^) 

C- z; j 2 (w L ) = er_ CT2 exp(-in7 2 • ip/Vk)(w L ) . (18) 



In order to find the wave function for general M, we 
note the following. The two bosons ip = (ipi,(f2) can 
be written in terms of a (hyper-)charge and spin boson 
(ip c and ip s respectively) by means of a simple rotation: 
ifi! = ^ip c + §</j s and tp 2 = -~<p c + ^ifis- The M- 
dependence is then brought in via a rescaling of the scale 
associated with the charge boson tp c . The particle and 
quasihole operators for general M become in terms of 
these bosons 



J 



k,M 

NASS.qh 



= lim (z oc ) a (exp 



B ai = V>i exp(-5=(V2fcM + 3<^ c + </? s ))(z T ) (19) 
v2k 



B 



fa exp(^(V2fcM + 3<^ c - (fi s ))(z L ) (20) 



1 



V2fcM + 3 
1 



Vc + ^ s )){w r ) (21) 
(p c -<p.))(v^) , (22) 



where we have written tp a 



ipt, -0-Q2 = 1P2 , cr roi = er T 
and <J_ ro2 = 17 1 for simplicity. The most basic spin fields 
at I transform as a doublet of the SU(2) subalgebra we 
identify with the spin of the particles. Also, the hyper- 
charge of the quasihole operators has the same sign as 
that of the particle operators. This implies that these 
are indeed quasiholes, as in earlier cases HQ]; the wave 
functions are nonsingular as any particle coordinate Zi 
approaches any quasihole coordinate Wj. 

Note that when these operators are used in the CFT 
correlator (together with a suitably chosen background 
charge), the extra Laughlin factor is automatically gen- 
erated. The correlator for the component of the wave 
function with iV-fj spin- up and down particles and n^j 
spin-up and down quasiholes is given by 



[V2kM + 3{N ] +N l ) + 



ip c + [iV T - Ni + ?i T - n^ifis \ (zoo) 



/2k I/' ■ - > ■ *' ■ ^/2kM + 3 

x C mi (w\) . . . C4 (w^ )C'_^ (w{)... C'_^ (w^ )B ai (4) . . . B ai (z^ )B_ a2 {z{) . . . B_ a2 )) 



(23) 

The value of a will be given momentarily. In the wave function (|23|) we inserted the most general background charge 
required for neutrality in the Cartan subalgebra of SU(3). However, the background charge can involve only the 
charge boson tp c , which corresponds to the spin-independent background magnetic field in the QH problem. Thus we 
find the condition 



(24) 



which is part of the requirement of SU(2) symmetry for the correlator. The correlator is a spin-singlct, which means 
that the wave function for the particles is a nonzero spin state of the particles, with spin determined by the quasiholes. 
Effectively, the quasiholes carry spin 1/2 which is added to the spin-singlet ground state. Note that a quasihole labeled 
up carries a spin down from the latter point of view, by Eq. (|24|), just as it carries negative charge (hence the term 
quasihole), since there is a deficiency of particles at its location. For N = iVf + iVj_ sufficiently large, N > n in fact, 
the spin (up or down) for each quasihole can be chosen freely, as we will see in some examples. Using condition (|24|), 
we can calculate that a must be 



a= 2W + 3 
2k 



2kM + 3 



(25) 



where n = n-\ in order that the limit z^ — > 00 exists and is nonzero. 

By working out the vertex-operator part, we arrive at the following form 



7 
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*§* A \*h »})] 1/k - W)* IN - 



n( 



w, — 10, 



)2W+3(f+ M )]J( W T 



y;* J 2feM+3 V fc 



r(f+M) 



(26) 



Note that the correlator is non-zero only if the 
parafermion and spin fields can be fused to yield the iden- 
tity operator. 

The number of magnetic flux N^, seen by any particle 
is found to be 

2fcM + 3 1 
^0= 2fc N+—n-(M + 2), (27) 



where we used Eq. (|24|). Since must be an integer 
(so that the wave function is a polynomial in the z^s), 
this gives another condition, that (3N + n) /2 [which is an 
integer by Eq. (^Jj)] must be divisible by fc. [For the RR 
states, there is an analogous condition, 2N + n must be 
divisible by fc. For fc even, this means that n is even, as 
in the fc = 2 case (the MR state) ||. In Ref. @, only the 
case n and N both divisible by fc was considered.] From 
Eq. we can deduce that the quasiparticle charge is 
l/(2fcM + 3). This corresponds to a fractional flux, l/2fc 
of the usual flux quantum. In effect, the flux quantum 
has been reduced by 1/k by the formation of clusters, as 
in paired states and in spin-polarized RR states j?J . The 
factor of 1/2 is present already in the Halperin k = 1 
case. So if fc is not divisible by 3, the quasihole charge is 
l/Q (Q is the denominator of the filling factor, defined 
in Sec. |n]), as in many other cases, but if fc is divisible by 
3, the quasihole charge is 1/3Q. This is similar to what 
happens in the MR and RR states, where the quasihole 
charge is further fractionalized (smaller than l/Q) when 
fc is divisible by 2 pf. 

The conditions (|2J) and that be integer are clearly 
necessary, but in fact are also sufficient, to ensure that the 
quasihole wave functions are nonzero polynomials in the 
ZiS, except in the special case n = 1 where the function 
vanishes. To see this, one must examine the fusion rules 
for the parafermion system, and check that the fields can 
be fused to the identity operator under the stated condi- 
tions. This will be considered in the next Section. 

For completeness, we give the conformal dimensions 
of the particle and quasihole operators B a (z) and 
C m (w). To obtain these, we need the dimensions of the 
parafermionic and spin fields, which are Aj, = 1 



A„ = 



k-l 
fc(fe+3) 



also 



A part 



respectively [|18| 
M±l and A, 



2 "q h 2(fc+3)(2fcM+3) ' 

We can show that the quasihole states we have ob- 



Using these, we find (see 

(5fc-l)A/+8 



tained are zero-energy eigenstates for the fc + 1-body 
Hamiltonian above, as follows. We again concentrate on 
the case M = 0. The argument using the ope's of the 
currents B^i(z) again applies M, as long as the fc + 1 
z^s are brought to a point that does not coincide with 
a quasihole coordinate Wj. To complete the argument, 
we must also consider the case where the latter condition 
does not hold. There are two ways to do this. One is 
to note first that, as a function of the Zj's for fixed w/s, 
the correlator is a polynomial, as it must be to be a valid 
QHE wave function. (It is not a polynomial in the quasi- 
hole coordinates Wj). This is because we chose to exam- 
ine quasiholes rather than the opposite charge objects. 
Then the fact that it vanishes when fc + 1 Zi's coincide 
away from a quasihole coordinate Wj also implies that 
it vanishes when they are at a Wj, by continuity, which 
holds because the function is a polynomial in the z^s. 

A second argument is also instructive. We may gener- 
alize the argument using the current algebra null vectors 
to directly address the limit of B's approaching a wj. 
There is a generalization of the central equation for this 
case, 



( J B T ,_ 1 ) fe |T) = o. 



(28) 



Here the state | f) = C Wl (0)|0) is the state in radial 
quantization corresponding to the quasihole operator at 
0. There are similar equations, with successively lower 
exponents, for the higher-order quasiholes (with multi- 
ples of the charge of the basic one) obtained by succes- 
sively fusing quasiholes together. These are the null vec- 
tor equations for the highest weights in distinct represen- 
tations of the affine Lie algebra (or for the distinct pri- 
mary fields) [^l],^) ■ We want to emphasize that the equa- 
tion says that for certain choices of the spins, the wave 
function vanishes when only fc particles come to the same 
point (or fewer for the higher-order quasiholes). This is a 
generalization of the fact that the Laughlin quasihole is 
defined as the factor ^(zi — w) which vanishes when any 
one particle approaches w. (This generalization applies 
already to the spinless RR states.) It is also a general- 
ization of the Halperin case, where a quasihole is a factor 



rw 



, r 



■jt) which vanishes when a single up particle 



goes to w* , but not when a single down particle does. 
When fc basic quasiholes are fused at w (taking the lead- 
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ing term at each fusion), the null vector equations state 
that the function does vanish when a single particle of 
appropriate spin approaches w, so we have a Laughlin- 
(or rather Halperin-) type quasihole in that case, as was 
already known for the MR (k = 2 RR) state, for example. 
Note that the above spin k/2 state is the highest weight 
in a multiplet, so there is a set of k + 1 such null vectors 
in total. This does not include all possible spin choices, 
as we have pointed out. To complete the proof that the 
function vanishes when any k + 1 ZiS (i.e. k + 1 particles 
of any spin) come to Wj, we must show that, for all m, 

(i? L _ 1 )™(i? T ,_ 1 ) fe + 1 - m |T>=o 

(B lrl r(5 trl )W- m |I} = 0. (29) 

This can be done by an elementary argument, applying 
the SU(2) lowering operator to Eq. (p8|), then another 
using the same equation, and then lowering fur- 
ther, and so on. 



VI. FUSION RULES 



From now on, we focus mainly on the case k — 2, which 
is a spin-singlet analogue of the MR state. The numeri- 
cal studies reported in Section [x| were all performed for 
this special case. Analytical results for k > 2 will be 
presented elsewhere p7f . 

As pointed out in the introduction, the non-trivial fu- 
sion rules play a crucial role in the ground state degenera- 
cies. In fact, the correlator in Eq. ( |26|) does not represent 
a single wave function, because in general there is more 
than one way in which the spin helds and parafermion 
fields can be fused to the identity. To show how this 
works, we will give the fusion rules, and explain that 
they can be written in terms of a Bratteli diagram. 
By using the correspondence between the fields of the 
parafermion theory and fields of the corresponding Wess- 
Zumino-Witten models (sec flq]), one finds the fusion 
rules listed in Table |. 

An examination of the fusion rules shows that there 
are different cases, according to the parity, even or odd, 
of raj, n j, iVf, iVj_. In the case where all four numbers are 
even, the spin fields and the parafermions can be fused 
to the identity separately. In the case where all four are 
odd, the spin fields and the parafermions can be fused to 
012 separately, and these two 0i2 's can then be fused to 
the identity. Because the quasiholes only involve the cr-p ^ 
fields, we in fact only need the first two columns of Table 
H With this restriction, the fusion rules can be written 
in terms of a Bratteli diagram, see Figure 




FIG. 1. Bratteli diagram for the spin fields of 
SU(3) 2 /[U(1)] 2 . 

Each arrow stands for either a cr-j- or field. The 
arrow starts at a certain field which can only be one of 
the fields on the left of the diagram at the same height. 
This last field is fused with the one corresponding to the 
arrow, while the arrow points at a field present in this 
fusion. As an example, the arrows starting at the * are 
encoding the fusion rules p x cr^n\ = "02(1) + a Ul) an d 
03 x (r«h = ^1(2) + o'xd)- One checks that the diagram 
is in accordance with the first two columns of Table Q. 
The symbol 1 at the right hand side of Figure |l| indicates 
that in the end we have fused the fields to the identity. 
This is possible only when nj and nj, are both even; in 
the case where both numbers are odd, 0i2 is obtained 
at that position in the diagram. In the remaining two 
cases, where n = n-\ + n± is odd, we can draw a similar 
diagram with the last point at the top, representing tpi 
or ip2 (except when n — I, to which we return in a mo- 
ment). In these cases, N = Nf + N± must also be odd, 
in order that the fusion of the parafermions 0i and 02 
for the particles can produce the appropriate field which 
can fuse with the result of the spin field fusions to finally 
give the identity. In the case n = 1, it is not possible 
to fuse the spin fields to obtain a parafermion, and the 
correlator vanishes. 

For the counting formula we need to know the number 
of ways in which a given number of spin fields can be 
fused to give a field in the parafermion sector, that is 1, 
"01 j 02, or 0i2. This number equals the number of paths 
(of length the number of spin fields) on the Bratteli dia- 
gram leading to the corresponding point on the diagram. 
One finds that when the numbers of spin-up and down 
quasiholes are 7ij and n±, respectively, this number, the 
number of fusion channels, is d„ Tjrii = f n -2, where the 
Fibonacci number is defined by f m = f m —x + fm-2, with 
the initial conditions /q = 1 and fi = 1. This is valid for 
11 > 2 in all four cases of rip, n± even or odd, while for 
n = l, d n ^^ ni = 0, and for n — 0, c?o,o = 1. The result is 
obtained by examining the Bratelli diagram and seeing 
that the number of paths obeys the recurrence relation 
that defines the Fibonacci number. We note that this 
result, the Fibonacci number f n -2, is the same as for n 
quasiholes in the k = 3 RR states 0. This is a mani- 
festation of level-rank duality, here between SU(3)2 and 
SU(2) 3 . 

The final fusion of the spin fields with the parafermions 
from the particles must produce the identity, in order 
that the correlator be nonzero. When d nr , ni is nonzero, 
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necessary and sufficient conditions for this are N-\ + n-f = 
N 1 + ji 1 (mod 2) and A T -n l = (mod 2). The first con- 
dition is the mod 2 version of the condition (^4j), while 
the second is equivalent to the condition (3N + n) /2 = 
(mod 2), on using Eq. (p4|). Note that, as can be inferred 
from the fusion rules for the tp fields (see Table |), the 
fusion of the i/j fields from the particle operators does 
not increase the degeneracy. The Abelian properties of 
the ip fields correspond to the Abelian statistics (Bose or 
Fermi) of the underlying particles. 

These results give the degeneracy d ni , ni of quasiholc 
states for fixed positions and spins of the quasiholes, 
which is the basis for non-Abelian statistics properties. 
We see that the result does not depend on how many 
quasiholes are spin up or spin down, and the sum over 
all choices of spins gives a further factor of 2™ spin degen- 
eracy for sufficiently large N, when only the positions are 
fixed. In the following Sections, we examine the total de- 
generacies of quasihole states when both their positions 
and spins are unrestricted. These are more suitable for 
numerical checks, and are finite numbers on the sphere 
(for a disk in the plane, they are infinite, and contain 
information about edge excitations as well as bulk quasi- 
holes). As in cases studied earlier 0,0, the total de- 
generacies are not just the numbers found above times 
a factor for the spatial degeneracy contribution, but in- 
volve partitioning the Fibonacci numbers above into a 
sum of positive integers. We also note here that when 
a generic Hamiltonian has a ground state in the NASS 
phase, the degeneracies will not be exact, but will be as 
given in this Section when all quasiholes are asymptoti- 
cally far separated. This will not be considered further 
in this paper. 



VII. SPATIAL DEGENERACIES 



Techniques for calculating degeneracies for a spheri- 
cal geometry are described in full detail in (l2). On the 
sphere, the relation between the number of particles and 
the number of flux quanta for the ground state is given by 
N$ — v~ 1 N—S. By increasing the number of flux quanta 
at fixed N, quasiholes are created. Moreover, when flux 
have been added, there may be zero-energy states with 
N not satisfying the conditions required in the ground 
state, for example that N be even in the MR state [fL2[ . 
In general, we would define the number of flux "added" 
as AN,/, — Nfj, — v~ 1 N + S. This is defined relative to 
a ground state at the same number of particles, even 
though such a zero-energy state for N not divisible by 
2k (or k for the RR states) would require a non-integer 
number of flux, and does not exist. Consequently, while 
our is always an integer, as discussed in Sec. [v|, AA^ 
does not have to be an integer. The number of quasiholes 
n can then be defined as n — kAN,p (for the RR states), 



or n — 2kAN c / 1 for the NASS states considered here, in 
agreement with our formula for N$ in Eq. (|27]) above. 

To explain the spatial degeneracies, we use the Laugh- 
lin case as an example, and follow the discussion of Jl2| . 
The Laughlin wave function for N particles in the pres- 
ence of n quasiholes can be written as fll| 



,w n ) 



(30) 



i<3 



t . k 



For this state, n = AiV^ (adding one quantum of flux 
creates one quasihole), so we have N$ = M(N — 1) + n, 
where we used that S = M for the Laughlin state. To 
continue, we calculate the degeneracy due to the presence 
of the quasiholes, by expanding the factor Y[ t k( z i ~ w k) 
in sums of products of the elementary symmetric poly- 
nomials 



E 



W^Wiz ...Wi„ 



(31) 



H<*2<"'<«n 



Viewing the coordinates w; as coordinates of bosons, we 
find that n bosons are to be placed in N + 1 orbitals. 
The dimension of the space of available states (linearly- 
independent wave functions) equals the number of ways 
in which one can put n bosons in ./V + 1 orbitals, which 
is 



N + n 
n 



(32) 



This is the spatial degeneracy we are after, although for 
the simple case of the Laughlin states. 

The situation for the MR state is discussed in detail in 
fll2| |. For the MR state, there is an additional complica- 
tion, namely quasihole states in which there are unpaired 
fermions are possible; this is the origin of the degenera- 
cies for fixed quasihole positions, already discussed. We 
will denote the number of unpaired particles by F, with 
the requirement that N — F be even, so that the num- 
ber of unbroken pairs (N — F)/2 is an integer. For n 
sufficiently large, N need not be even in the zero-energy 
states. It was found that the spatial degeneracy depends 
on the number of unbroken pairsx; in fact, for the MR 
state, it was given by |Q 



N-F 



(33) 



For the clustered state of [Q (in which the particles 
form clusters of order k rather than pairs), the spatial 
degeneracy is given similarly by j7 13 



N-F 



(34) 



where N — F must be divisible by k (there appears to 
be no known general analytic proof of this formula for 
k > 2). 
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Based on these earlier results, and because the NASS 
wave functions involve clusters of up particles and sep- 
arately of downs, we expect that the spatial degeneracy 
for the NASS states is just a product of two binomial 
coefficients, involving the -Fi, F2 "unclustered" particles 
(or parafermions) for spins up and down, respectively: 



vanishes if the conditions that Fi 



JV T -Fi 



N1-F2 



(35) 



Again, iVj— F\ and N^ — F 2 must be divisible by k. Notice 
that these numbers never depend on M. 

The case k = 1, F± = F2 = gives the spatial de- 
generacies for the general Halperin (to', to', to) Abelian 
states (as was mentioned briefly in Ref. for a partic- 
ular case). For the spin-singlet cases (m + l,m + l,m) 
(to = M + 1) of interest here, the conditions that cor- 
respond to zero-energy states with only the particle and 
quasihole numbers, N and n, fixed are that + = 
Ni + rij. For N sufficiently large, these allow any choice 
of spin, up or down, for each quasihole, giving a factor 
of 2™ spin degeneracy for fixed positions (this holds for 
all k). Such degeneracy does not contribute to the de- 
generacy on which non-Abelian statistics depends, and 
the statistics is Abelian in the present k = 1 case (as of 
course was expected). The full degeneracy in this case 
k = 1 is obtained by summing ( |35| ) over jVf, JVj_, n-f, n^, 
subject to the conditions just stated, with F\ = F 2 = 0. 
These imply that the summation is over only the possible 
values of S z = (N^ — JVj,)/2, and the result is 



N ■ 



(36) 



Note that this number includes the spin degeneracy. If 
we take the ratio to the number 



{N + n)/2 



(37) 



of quasihole states with, say, n± — 0, the result tends 
to 2™ as N — > 00, which is again the spin degeneracy 
for fixed positions. In the remainder of this paper, we 
concentrate exclusively on k = 2. 

With the orbital degeneracies in hand, we need to know 
how to break up the degeneracies stemming from the fu- 
sion rules of the spin fields. So in fact we need to know 
the number of unpaired particles of either spin, Fi, F 2l 
for each possible path on the Bratteli diagram. The next 
Section will treat this problem, by partitioning the num- 
bers d n] . iTll in the following way 



E 

■2*1 



F\ F 2 



(38) 



The symbol {jp p l }i is interpreted as the number of zero- 
energy states containing n-\, nj quasiholes at fixed posi- 
tions and F\, F 2 unpaired parafermions. The symbol 



F 2 — rij be even 



are not satisfied; these conditions are equivalent to the 
two mod 2 conditions discussed at the end of Sec. VI. 



The subscript 2 indicates that we are dealing with the 
case k = 2. 



VIII. COUNTING SU(3) 2 /[U(1)] 2 PARAFERMION 
STATES 



It was explained in jl2| that the state counting for the 
MR state involves the systematics of Majorana fermions, 
which act as BCS quasiparticles (unpaired fermions) oc- 
cupying zero-energy states |20| . For the more general 
RR (spin-polarized) states with order k clustering [Q, 
the Majorana fermion is replaced by an 811(2)^/11(1) 
parafermion 0JH]- 

We recall that for the NASS state at level k = 2, 
the role of the Majorana fermion is taken over by the 
parafermions that are associated to SU(3)2/[U(1)] 2 . The 
(spin-up or down) quasiholes correspond to (two differ- 
ent) spin fields <7f.| of this parafermion theory (see Sec- 
tion |V|). In Section |\T|, 



we calculated the number of 
different quantum states (conformal blocks for the corre- 
lators) that can result from introducing n-f cr| and nj 
spin fields (and also varying the number of particles N-\ , 



Ni) with the result d n 



f, 



-2- The degeneracy 



results from the presence of varying numbers of particles 
that are not members of clusters, which in the correla- 
tors can be identified with the fundamental parafermions 
ipi, ip 2 of the parafermion theory. F\ (F 2 ) is the num- 
ber of tpi (1P2) excitations. These numbers are subject to 
the conditions that F\ = uj, (mod 2), F 2 = »f (mod 2), 
otherwise the number of zero-energy states is zero; these 
conditions come from those discussed in Sec. VI. In the 



previous Section, we found that the orbital degeneracy 
depends on the numbers Fi and F 2 , see Eq. (|35|). Wc 
now turn to the calculation of the symbols 



Fi F2 



(39) 



which partition the degeneracies d nr ^ ni in the correct 
way. We will also keep track of the angular momentum 
(L) multiplet structure associated with these parafermion 
states. To do this, we have to go through a series of steps. 

First of all, we consider the infinite (chiral) character 
corresponding to the full parafermionic CFT (see [p5[) 



ch(xi,x 2 ;q) = 



Fi,F 2 



q { F ? +F *- FlF2)/2 x F lx F 2 
(q)f 1 {q)F 2 1 2 



(40) 



where (q) a = \~[j=i(^ ~ 9*0 f° r integer a. Here Fi and F2 
are unrestricted non- negative integers, and x±, x 2 , q are 
indeterminates. 
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What is needed for our purposes here is the trun- 
cation of this expression to (a sum of) polynomials 
Y n f,m (xi,x 2 ; q) that describes the states that occur when 
n|, ni spin fields (quasiholes) are present. The approach 
is described in Refs. ]l4],^6|], see also Ref. Jl3[, and details 
for the present case will be given in Ref. ]27|]. We find 
that these polynomials satisfy the following recursion re- 
lations 

Y( nr>ni ) = i / („ T -2,n i ) + X\q 2 i / (n T -2,n i +l) , 

"j-i 

Y(nt,ni) — ^(n T ,n i -2) + X 2<1 2 Y(n r +l, ni -2) (41) 

with initial conditions 

Y(i,o) = ^(o,i) = > 

Y(Q,0) = ^(2,0) = ^(0,2) = 1 , 

F (M) = qix x x 2 . (42) 

Recursion relations similar to the above (but lacking the 
xi.2 dependence), have been considered in the mathe- 
matical literature on special polynomials associated to 
SU(3)2, see for instance |^| . The coefficient of x Fl x^ 2 



Y(8,0) 
Y(7,l) 
F (6,2) 

y (4,4) 

etc. 



in the polynomial Y"( nTin ,) is a polynomial in q with the 
sum of the coefficients equal to the symbol ([39]), that is 

Y (nuni) ( Xl ,x 2 ,l) = Y < lx 2 {f F 1 } ■ (43) 

F U F 2 ^ 1 2 ^ 2 

We notice that the recursion relations preserve the con- 
ditions on the parities of F\, F 2 (the exponents of x\, X2) 
that are part of the definition of the symbol (^). In the 
limit where (n-j-, n\) — > (00, 00), the sum of these polyno- 
mials over the four choices, n-f and nj, each either even or 
odd, approaches the expression ch(xi, x%; q) given above. 

The coefficient of xf 1 x F2 in the polynomial Y n ^ <n . has 
a special form, which allows us to extract information 
on the L quantum numbers of the parafermion states: 
after multiplying with a factor g - ( n T- F i+ n 4.- F i>)/ 4 ) W e rec- 
ognize a sum of terms of the form q z , which together 
form a collection of angular momentum [L) multiplcts 
with quantum numbers L z — l z fig] ]. 

To illustrate the above, we present the polynomials for 
the case of two added flux quanta, giving eight quasiholes. 
The polynomials are 



= 1 + (q 2 + q 3 + 2q 4 + q 5 + q 6 )xj + q S xf + (q 6 + q 7 + q 8 + q 9 + q W )x\xl , 

, 1 3 5 7 T .,9 11 13 15. o 19 r: q 

= (q 2 +q 2 +q 2 + q 2 )x 1 x 2 + (q 2 + 2q? + 2q 2 + 2q 2 + q 2 ) X fx 2 + q 2 x\x 6 2 , 

= 1 + (q 2 + i + q^xf + (q 2 + q 3 + 2q 4 + q 5 + q*)x\x\ + (q e + q 7 + q*)x\x 2 2 , 

= (g3 + 2q% + 2q% + q^)xix 2 + + q^ 2 + q^)x\x 2 + (q* + q^ + q^ + q^)x\x\ , 

= 1 + q 2 x\ + q 2 x 2 + {q 2 + 2q 3 + 3q 4 + 2q 5 + q 6 )x 2 x 2 + q S x\x\ , 

(44) 



From the polynomial ^(5,3) (as an example), we read off 
the following nonzero values of the symbols 



5 3 
1 1 

5 3 
3 1 

5 3 
3 3 



3 1 

6 (L = -,L=-) 
1 2 ' 2' 



3 (L = l) 



= 4 i L =i> 



(45) 



In fact, it is possible to write the polynomials Yr n . tn ,^ 
in a closed form p9[ , 

Y (mtni) ( Xl ,x 2 ;q) = Y, 'q {Fi+Fi - F ^ ),2 x F ^x^ 

Fi,F 2 



[ n r +F 2 - 




1" ni+Fi 1 


2 




2 


^1 




F 2 



where 
fined as 



, (46) 

is the g-deformed binomial (q-binomial), de- 
= ( q )lfq) — j) an d the prime on the sum 



denotes the restriction on F±, F 2 values. Using the prop- 
erty that in the limit q — > 1 the g-binomials become or- 
dinary binomials, we find the following explicit formula 
f° r {}2 (under the same conditions on Fi, F 2 , otherwise 
it vanishes): 



Fi F 2 



n^+F 2 
2 

Fi 



2 

F 2 



(47) 



Note that if we take the sum over all F\ and F 2 , we indeed 
find the correct value, namely a Fibonacci number 



, / n r +F 2 
\ I 2 



Fi ,F 2 



n t +Fi 
2 

F 2 



fn^+ni —2 



(48) 



While Eq. ( f47| ) gives us just the number, we can also 
obtain the angular momentum content easily. The bi- 
nomial ( °j \ is interpreted as the number of possible 

ways of putting / fermions in a boxes which have quan- 
tum numbers L z = —(a— l)/2, ... , (a— l)/2 assigned to 
them. In this way, an angular momentum multiple struc- 
ture is assigned to the binomials (see Q). The angular 
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momentum content of the symbols {}2 is obtained by 
adding the angular momenta associated to the binomials 
in the usual way. 



IX. FINAL COUNTING FORMULA 

We are now in the position to write down the formula 
for the total degeneracy of zero-energy quasihole states 
of the k — 2 non-Abelian spin-singlet states. Recall that 
there are two conditions on the numbers of quasiholcs 
(see Section |v|). The first condition is N^+n^ = N± 
which implies that the correlator is a spin-singlet, or that 
the wave functions have total spin determined by the 
spin- 1/2 quasiholes added. The other condition was that 
(37V + n)/2 be even, to ensure that is an integer, 
where N = iVj + iVj_, and n = nj + ni = 4AiV0, which 
relates the number of excess flux quanta and the number 
of quasiholes added. These imply that iVf — nj. = N± — n-f 
must be even. 

The result of the previous few Sections is now that the 
total number, summed over all spin components, of zero- 
energy states as a function of the number of particles and 
added flux quanta is 



':(N,ANt) = 

Nf-F! \ / N l -F 2 



n,+F 2 
2 



n l +F 1 
2 

F 2 



n j 



2 ' n i 



(49) 



where the prime on the sum indicates that it is restricted 
to values obeying all the conditions just mentioned, and 
to iVf — Fi and 2Vj — Fi even as discussed in Sec. VTJ. 
Note again that these conditions imply that n-j + F% and 
nj, + F\ are even. 

In addition, the orbital angular momentum decompo- 
sition of the states can be obtained, by combining the 
angular momenta found in the orbital and parafermion 
factors in the preceding two Sections. The spin quantum 
number of any given state is simply S z = (JVf — iVjJ/2 
and one readily recognizes the multiplet structure for the 
SU(2) spin symmetry. (We remark that the parafermion 
theory by itself does not have a proper SU(2) spin sym- 
metry.) 

In Table [n], we present counting results for N = 4, 8, 
12 and AiV^ = 1, 2, 3, 4. We specify the number of states 
as a function of the L and S quantum numbers. In Ta- 
ble III we give some results for N not a multiple of four. 
Notice that for n = 1, there are no zero-energy states, as 
expected from Sec. VI. The results listed in Table || and 



[II are for the cases we checked numerically, as we discuss 
in the next Section, and are in full agreement with those 
results. 



X. NUMERICAL METHODS AND RESULTS 

We next turn to some numerical studies of the NASS 
states. We consider only cases where the particles are 
fcrmions, to represent electrons. We have studied the 
k = 2,M = l(u= 4/7) case in both the toroidal (PBC) 
and spherical geometries. We first present the results for 
the sphere. As discussed before the flux-charge relation 
for this state is — 7N /4 — 3. The number of single- 
particle orbitals (the lowest LL degeneracy) is N$ + 1 . In 
order to make contact with the results on more conven- 
tional geometries the radius R of the sphere has to be 
chosen so that the number of flux is N<j, = 2R 2 (where 
the magnetic field strength B is fixed, such that the mag- 
netic length is 1 in our units), so R = y/ N^/2 p6| . The 
filling factor is v = N/N^ = 2nn, where n = N/(4ttR 2 ) 
is the particle number density. 

For numerical purposes, it is best to re-express the in- 
teraction Hamiltonian in terms of projection operators 
onto different values of the total angular momentum for 
different groups of particles |l6|. For the M = 1, k = 2 
case of the NASS states, the required Hamiltonian can 
be written as 



H = U Puk$N+/2 -3, 3/2) + V' 

i<j<k 



i<3 



(50) 

with U, V > 0. Here P ljk {L,S) (P i3 (L, S)) are projec- 
tion operators for the three (resp., two) particles spec- 
ified onto the given values of total angular momentum 
L and spin S for the three (resp., two) particles. Each 
projection is normalized to P 2 = P. To see that this is 
the required Hamiltonian, that corresponds to the short 
range 5-function interaction for M — 0, and gives the 
same numbers of zero-energy states found above, note 
the following. First, the maximal angular momentum for 
several particles corresponds to the closest approach of 
those particles Jlq] , In particular, the two-body term is 
a contact interaction, and V — Vq, the zeroth Haldane 
pseudo-potential fig] . The two-body term implies that 
any zero-energy states must have no component with to- 
tal angular momentum and total spin zero, which, 
since we are dealing with spin 1/2 fermions, means the 
wave function must vanish when any two particles coin- 
cide. The wave function must therefore contain a factor 
multiplication by this factor defines a one-one map- 
ping of the full space of states of spin 1/2 bosons in the 
lowest LL, with reduced by N — 1, onto the sub- 
space of states of the fermions that is annihilated by the 
two-body term in H. Under this mapping, the three- 
body Hamiltonian for the M = case corresponds to 
the three-body term in H, and selects the correspond- 
ing states as zero-energy states. In particular, the total 
spin of the three bosons when they coincide (and hence 
of the fermions) must be 3/2. Hence the zero-energy 
eigenstates of the present Hamiltonian are given by the 
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results derived earlier. Note also that H can be rewrit- 
ten in terms of ^-functions and their derivatives. The 
zero-energy eigenstates of this Hamiltonian were found 



for various N and values, and analyzed in terms of L 
and S. The results are shown in Tables || and III, and 
agree with the counting formulas presented above. 
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FIG. 2. The spectrum of the NASS model ground state for N = 8 and 4/7 filling. The last panel shows all S values combined. 
The insets are the low lying levels. 



Next we discuss the full spectra of the Hamiltonian. 
In Fig. ^ we show the excitation spectrum of an N = 8, 
N<f, = 11 system classified by the total spin S = to 
5 = 4. Whenever necessary we have shown the low-lying 
spectrum in an inset. The frame in the lower right hand 
corner shows the entire spectrum irrespective of the to- 
tal spin quantum number of the state. The choice of U 
and V is immaterial to the ground state, which is always 
the unique zero-energy eigenstate of H. Obviously, the 
excitation spectrum will depend on the choice of these 
coefficients. In producing Fig. [| we chose U = V = 1 so 
that H is a sum of projection operators. There appears 
to be a well-defined gap, suggesting that the system is 
incompressible (for spin as well as charge) in the ther- 
modynamic limit, as assumed in the preceding analysis. 



In this connection, we may point out that, as well as 
the quantized Hall conductivity for charge, our system 
then has a quantized Hall conductivity for spin, given 
by fc/47r in natural (h = 1) units, which is associated 
with the SU(2)fc subalgebra of the chiral algebra (see 
e.g. Ref. 1 20 1 and references therein). Collective modes 
with S = [L = 2,3,4) and S = 1 (L = 1,2,3) be- 
low the continuum can be tentatively identified in the 
spectra (see insets in Fig. ||). That is, these may be 
finite-size dispersion curves of single neutral excitations 
in plane-wave (spherical harmonics on the sphere) wave 
functions, which would be charge and spin modes, respec- 
tively. We shall not address the precise nature of these 
neutral modes here. 
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ground state of the Coulomb potential. However, our 
numerical data for N — 8 shows a much stronger cusp 
at the N,), of the NASS state than at the of the state 



FIG. 3. The spin up-up and spin up-down pair correlation 
functions, together with their sum (solid line) and difference 
(dashed line), versus the chord distance, calculated in the 
ground state for TV = 8, iV^ = 11 (y = 4/7). 

In Fig. H we show the various pair correlation func- 
tions of interest, g^(r), gn(r), as well as <7Totai = 
<7ff (r) + 9]i(r) and 5tl ( r ) — $tt( r )- The widely-different 
correlations between like and opposite spins is no doubt 
magnified by finite-size effects. 

The 8-electron system is the first non-trivial size and 
is probably too small for any meaningful comparisons 
of the overlap with the 2-body Coulomb potential prob- 
lem. Nonetheless we found a nontrivial overlap-squared 
(about 55%) with the ground state of the Coulomb poten- 
tial for particles in the lowest LL (again, with no Zeeman 
term), at the same N, N^. By modifying the value of the 
lowest pseudo-potential for the Coulomb interaction this 
overlap-squared can be improved to 93% (and probably 
beyond) without any intervening phase transition (an en- 
ergy gap with the ground state is maintained at all times 
while the pseudo-potentials are varied). However, in the 
lowest LL we do not expect to produce a better trial wave 
function than one constructed by the composite fermion 
(CF) method fl3(i|| , in which two flux quanta per particle 
are attached (in the opposite direction to the background 
magnetic field) and the resulting CFs fill completely the 
first two LLs of CFs (with both spins). By construction 
this is a uniform (L = 0) spin-singlet state. We have 
not constructed this state, as it occurs at a different flux 
for a given N (N^, = 7N / 4) , making a direct comparison 
with our NASS state difficult. We note, however, that for 
N = 8 the CF state corresponds to a spin-singlet Fermi- 
liquid-like state, as at v — 1/2. That is because the net 
effective field of the CFs is zero for this size (states that 
lie in sequences for different filling factors can coincide 
at finite size on the sphere, because the shifts S may be 
different — see e.g. [[H)). We expect that, as usual, this 
CF state will have a very large overlap with the exact 



obtained by the hierarchy/CF construction, where in fact 
no cusp can be discerned. See Fig. ^ below. 
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FIG. 4. The ground state energy of the pure Coulomb po- 
tential in the lowest Landau level versus N^, at N = 8. The 
numbers of flux N$ for the NASS and the spin-singlet CF 
states are marked. 



We have also studied the N — 8 size on the torus. Un- 
fortunately, as on the sphere this size is too small for any 
meaningful comparison (e.g. there are only four distinct 
many-body k vectors for this size; one is at the zone cen- 
ter, the other three at the zone boundary). We would just 
like to point out that, for the analog of H on the torus, 
the degeneracy for the 4/7 state is 2 (excluding the 7- fold 
center of mass degeneracy) , in agreement with the count 
in Sec. [b], since the number of particles is divisible by 
4. These are two k = states. For the pure Coulomb 
potential in the lowest LL the state at 4/7 is in fact com- 
pressible: The total spin is S = 1 and its k vector varies 
with the geometry of the PBC unit cell. However, one 
obtains an incompressible state by increasing Vq or V\ 
and we obtained overlap-squared as large as 50% when 
we compared the lowest two states (which happen to be 
both S = 0, k = states) to the model NASS states. 
Note that the shift is zero on the torus, and there can 
be no interference from v — 1/2 here. We have not per- 
formed any further or systematic studies of such issues 
because, as in the case of the sphere, we suspect that the 
CF-based state will be closer to that of the Coulomb po- 
tential. That is, we expect that the system with Coulomb 
interaction in the lowest LL at v = 4/7 is in fact in the 
hierarchy/composite-fermion phase (whether spin-singlet 
or not), not the NASS phase considered in this paper. We 
will return to more comprehensive studies of larger sizes 
in the future. 
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TABLE II. Counting results for the NASS states at k = 2. 
A^ is the number of electrons; AA^ is the number of excess 
flux quanta. The results are given as a function of the L (an- 
gular momentum) and S (total spin) quantum numbers. The 
total number of states is also indicated. 
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TABLE III. Counting results for the NASS states at k = 2 
with fractional AAV (symbols as in O) . 
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